The phonon dispersion of graphite by inelastic x-ray scattering 
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We present the full in-plane phonon dispersion of graphite obtained from inelastic x-ray scattering, including 
the optical and acoustic branches, as well as the mid-frequency range between the K and M points in the 
Brillouin zone, where experimental data have been unavailable so far. The existence of a Kohn anomaly at the 
K point is further supported. We fit a fifth-nearest neighbour force-constants model to the experimental data, 
making improved force-constants calculations of the phonon dispersion in both graphite and carbon nanotubes 
available. 



PACS numbers: 63.20.Dj, 63.70.+h, 78.70.Ck 

I. INTRODUCTION 

Research on carbon nanotubes and the recent availability 
of single graphite sheets (graphene) has revived the inter- 
est in the fundamental physical properties of graphite during 
the last years. Carbon nanotubes can be regarded as one or 
more rolled-up graphite sheets, and many physical properties 
of carbon nanotubes are closely related to those of graphite.-^ 

The fundamental characteristics of a crystalline material 
comprise its phonon spectrum, from which one can derive sev- 
eral other physical properties such as sound velocity, thermal 
conductivity, or heat capacity. Furthermore, phonons play an 
important role in excited-state dynamics and electrical trans- 
port properties. Optical or electronic excitations can decay 
into vibrational excitations or can be scattered by phonons 
into different states. For example, in carbon nanotubes the 
high-bias electrical transport is assumed to be limited by scat- 
tering of the carriers by optical phonons corresponding to the 
graphite K point. ^'^ 

The phonon dispersion of graphite has not been com- 
pletely resolved by experiment, mostly due to the lack of 
large enough samples of crystalline quality. It has been partly 
measured by inelastic neutron scattering (INS), electron- 
energy loss spectroscopy (EELS), and inelastic x-ray scatter- 
ing (IXS)i2i2i24i2iii Most experiments so far have determined 
the dispersion along the F — i^T and the F — M directions 
in the graphite Brillouin zone (see Fig. [T] for a definition of 
the Brillouin zone). The recent measurement of the optical 
branches along Wxq K — M direction by IXS pointed to the 
existence of a Kohn anomaly for the highest phonon branch 
at the K point J^^i^ Although this result resolved previous dis- 
crepancies between different lattice dynamics models, there 
are still open questions regarding the shape of lower-lying 
phonon branches. In particular, differences appear between 
force-constants and density-functional theory (DFT) calcu- 
lations, where experimental data are still unavailable. This 
concerns, e.g., the crossing between the acoustic and optical 
bands near the M-point or the energy of the transverse acous- 
tic mode at the K point. For carbon nanotubes, the experimen- 



tal determination of the phonon dispersion throughout the en- 
tire Brillouin zone would require monocrystalline samples of 
a minimum size, which have been unavailable so far. There- 
fore, the closest approximation to the experimental phonon 
dispersion of carbon nanotubes is currently the phonon dis- 
persion of graphite. 

Here we present the phonon dispersion of graphite in 
all three high- symmetry directions in the basal plane de- 
termined by inelastic x-ray scattering. In particular, the 
phonon branches between the K and M point and the acoustic 
branches in all high- symmetry directions are obtained, giving 
both the optical and acoustic phonons from one experimental 
technique. We fit our data by a set of force constants, includ- 
ing fifth-nearest neighbours of carbon atoms. The fitted force 
constants can be used to deduce the corresponding force con- 
stants for carbon nanotubes. 

This paper is organized as follows. In the next section we 
briefly describe the experimental details of the IXS experi- 
ments. We give an introduction to the phonon dispersion of 
graphite and present the experimental data in Sect. [Till In 
Sect. HV] we apply a fifth-nearest neighbour force constants 
fit to the experimental data and provide the in-plane and out- 
of-plane force constants. 



II. EXPERIMENTAL SETUP 

The inelastic x-ray experiments were performed at beam 
line ID28 at the European Synchrotron Radiation Facility 
(ESRF). For a review of IXS the reader is referred to Refs.[l3l 
andH The energy of the incident radiation of 17794 eV 
was selected by the (999) Bragg reflection of a silicon crys- 
tal. The scattered photons were analyzed by five analyzers 
operating in the same reflection order. The total energy reso- 
lution in this configuration is 3.0 me V.— The x-ray beam was 
focused to 250 x 60 /im^, selecting a single microcrystal in 
a naturally grown graphite flake. The typical size of a sin- 
gle grain was about 800 jim in lateral direction and 100/im 
along the c-axis. By x-ray diffraction we obtained the lattice 
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parameters a = 2.463 A and c = 6.712 A, in excellent agree- 
ment with previous neutron diffraction data (a = 2.464 A, 
c = 6.711 A).i^ 

Inelastic scattering spectra were recorded by varying the 
temperature difference between the monochromator and the 
analyzer silicon crystal. To minimize the effects of tempera- 
ture drifts that could result in an energy offset, we performed 
systematic Stokes-anti- Stokes scans between the measure- 
ments. In our setup the c-axis of graphite and the scattering 
plane encompasses an angle of 90°, 30°, and 0°, depending on 
the phonon branch under consideration. The scattering geom- 
etry was chosen according to the selection rules, see Ref . [l6[ 



III. EXPERIMENTAL RESULTS 

The unit cell of graphene contains two atoms, resulting 
in six phonon branches. The unit cell of graphite consists 
of four atoms, which leads to twelve phonon branches. The 
space group of graphite is P63/mmc (international notation). 
At the r point it possesses the factor group 6/mmm {Dq^ 
in Schonfliess notation). The optical zone-center modes of 
graphene are decomposed into T = B2g + E2g . In graphite, 
the optical zone center modes are decomposed into F = 
A2u + 2B2g + Elu + 2E2g jU^lMl The A2u and Eiu modes are 
IR active, the E2g modes Raman active. The B2g modes are 
optically inactive, but can be measured via INS or IXS. The 
three acoustic modes are decomposed into T = A2u + Eiu. 

Graphite is a highly anisotropic material: the nearest- 
neighbor distance between two atoms in the plane is a/ \/3 ^ 
1.42 A, while the inter-layer distance is c/2 ^ 3.35 A. The 
bonds between two carbon atoms in the plane are much 
stronger than the weak van-der-Waals interactions between 
the layers. Therefore, compared to graphene, one expects that 
the phonon modes of graphite correspond approximately to in- 
phase and out-of-phase vibrations of the two graphene planes. 
Most of the phonon branches in graphite are nearly doubly de- 
generate and almost the same as in graphene. ^^'^^ Only close 
to the r point, the acoustic modes of the single layer split 
in graphite into an acoustic mode (in-phase vibration of the 
graphene sheets) and an optical mode [out-of-phase vibration; 
in-plane: E2g at 5.2 meV (42cm~^); out-of-plane: B2g at 
15.7 meV (127 cm~^)]. For the optical modes of graphite, the 
difference between the in-phase and the out-of-phase vibra- 
tions is very small: at the F point the IR active Eiu mode 
is found at 196.9 meV (1588 cm~^), close to the Raman ac- 
tive E2g mode at 196.0meV (1581 cm~^). The same holds 
for the A2u mode at 107.5 meV (867 cm-^)i2iii and the B2g 
mode at 107.6 me V (868 cm~^). Therefore, in the following 
theoretical discussion, we will consider the phonons of a sin- 
gle graphene sheet. 

The six branches are divided into the out-of-plane acoustic 
mode ZA, the in-plane acoustic mode TA (sometimes called 
SH=shear), the longitudinal acoustic mode LA, the out-of- 
plane optical mode ZO, the in-plane optical mode TO (SH*), 
and the longitudinal optical mode LO. Four branches belong 
to modes where the atoms move in-plane with the graphene 
layer (TA, LA, TO, LO); two branches belong to transverse 




FIG. 1 : Brillouin zone of graphite. The Brillouin zone of graphene is 
the hexagon which lies in the plane with the points F, K and M. The 
distances between the high-symmetry points SLveF — K = 47r/3a, 
F - M = 27T/VSa, mi&K-M = 27r/3a. 



modes, where the atoms move out of the plane (ZA, ZO). 

In Fig. [2] we show our experimental data of the graphite 
phonon dispersion in the plane. The lines show the fifth- 
nearest neighbor force-constants fit described in Sect. [IVl 
The optical phonon frequencies near the F point agree well 
with previous experiments. We find the E2g LO mode at 
196.0meV (1581 cm-^), and the B2g mode at 107.6meV 
(868 cm~^). Regarding the overall shape of the phonon 
branches, our experiments confirm previous ab-initio DFT 
calculations, letting aside the special situation for the highest 
branch at the K point.Jii2i22i22i2^ 

As can be seen, the highest optical frequency does not ap- 
pear at the F-point. Instead, the phonon frequency first in- 
creases with larger wave vector and then decreases again. This 
effect, called overbending, has been observed in diamond as 
well.^"^ In graphite, it has been predicted to result from a Kohn 
anomaly, i.e., the frequency at the F point is lowered due 
to interaction of the phonon with the electronic system. 
Another Kohn anomaly in graphite can be found for the TO- 
derived phonon branch at the i^-point (fully symmetric A'^ 
{Ki) mode). We have gained additional data for the high- 
est optical phonons around the K point, confirming previous 
measurements of the frequency softening of the TO-derived 
branch near the K point. Again, we were not able to detect 
the A'l phonon directly at the K point. The strong electron- 
phonon interaction has been predicted to reduce the phonon 
lifetime which results in a line broadening. Probably the large 
line width makes it very difficult to detect the A'^ phonon at 
the K point experimentally. 

Considering the differences between previous theoretical 
models, we find the following results, see also Sect. [IVl Be- 
tween the F and M points, the ZO and TA modes do not cross 
within our experimental error of 3 meV. This is in contrast to 
previous empirical force-constants models and EELS data of 
Ref. 8. Which branch is higher directly at the M point cannot 
be uniquely distinguished from our data. However, in DFT 
phonon calculations the crossing is found between K and M 
but close to M (about 1/10 of the distance between K and M, 
see Ref. 23). The overall agreement with DFT calculations 
supports the crossing between K and M, see also Fig.O 

The TA branch along F — M shows a smaller increase com- 
pared to the electron energy loss spectroscopy (EELS) data in 
Ref. ^ In a recent EELS experiment on epitaxially grown thin 
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FIG. 2: (Color online) Phonon dispersion of graphite from inelastic x-ray scattering (symbols). Triangles are present data, circles are data 
already published in Ref. 11. The full square at the F point is INS data from Ref 7. Solid lines are the force-constants calculations from the 
5th-nearest neighbor fit discussed in Sect.|lVl the dashed line is a quadratic extrapolation of the data. The lines are denoted by their symmetry 
representation in space group notation. The relation between space group and molecular notation can be found in Table U 



We measured for the first time the ZA and TA mode be- 
tween the K and M points. Our results confirm predictions 
made by ab-initio calculations very well (Fig. [3]), and are also 
well reproduced by our force-constant fit (Fig. O. The trend 
of both branches crossing near the M point can be recognized. 

The two optical phonons at the M point derived from the 
LO and LA branches are very close in frequency (< 4 meV ), 
and we were not able to distinguish them clearly by symmetry. 
It appears, however, consistent with DFT calculations and the 
force constants fit in Sect. |IV] that the higher frequency has 
symmetry and the lower one . As a consequence, 
the LO- and LA-derived branches cannot cross between the 
K and M point. 

In Fig. H] we show the low frequency phonon range along 
the F — A direction, i.e., perpendicular to the in-plane direc- 
tion. For comparison, we also present the INS data on highly 
oriented pyrolytic graphite from Ref. 7. They are in excellent 
agreement. The high-frequency phonon range is expected to 
show almost no dispersion along the V — A direction. 

Due to experimental reasons we were not able to record data 
points from the ZO branch along the V — K — M direction. 
This branch has been measured in recent EELS experiments 
of Ref. [13 In general, the data of Ref. [Tol agree well with ours, 
but at the i^T-point the ZO and ZA branches in Ref. [13 show a 
relatively large splitting of ?^10meV. They cannot stem from 
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FIG. 3: (Color online) Our in-plane phonon dispersion of graphite 
together with a DFT calculation from ReffTTI (solid lines). Same 
symbols as in Fig. [21 were used. 



graphene sheets,-^ however, this branch could not be detected, 
as the shear modes in graphene are forbidden in EELS. These 
contrasting results suggest that the crystalline quality of Ref. 8 
was lower, softening the selection rules. This explains why 
some previous empirical models, relying on the then available 
EELS data, predicted a larger slope and consequently a cross- 
ing of the ZO and TA modes between the F and M points. 
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FIG. 4: (Color online) Phonon dispersion of graphite along the F — 
A direction. Open triangles are present IXS data, full squares are 
neutron scattering data from Ref.|7|. The lines are a guide to the eye. 



the degenerate i^e-phonon, but could possibly represent the 
out-of-phase modes of the graphite planes. On the other hand, 
in DFT calculations of graphite^^ this splitting seems much 
smaller than indicated by the EELS data. 

Regarding the phonon modes specific for graphite with 
more than one layer, we find the low-energy out-of-phase 
modes near the F point. These are indicated by the dashed 
line in Fig.[2l We measured two out-of-phase ZO' phonons in 
the F — M direction, with energies 13.6meV and 23.3 meV 
at 0.16 of the F - M distance and at 0.4 F - M, respectively. 
A quadratic extrapolation leads to a value of 12.5 meV at the 
F point, in agreement with 15.7 meV from neutron scattering 
data.i^ 

The optical-phonon frequencies at the high- symmetry 
points from our experiment are summarized in Figs. [6j U\ 
and [8] together with the displacement patterns obtained by the 
force-constants calculations in Sect. |lVl The acoustic phonon 
branches near the F point give information on the elasticity of 
graphite, which will be reported elsewhere.— 



IV. FORCE CONSTANTS CALCULATIONS 



FIG. 5: Coordinate system for the point group C2v of F — M and 
r — K drawn in real space, z points in the direction of the principal 
rotational axis. 



phonon branch is higher than the LA phonon at the M-point, 
and, as a result, the two branches do not cross between M 
and K; vice versa in most predictions by empirical force con- 
stants. 

Discrepancies with the experimental data are found for both 
models for the TO-derived branch at the K point, except for 
Ref. S In this context we want to emphasize the impor- 
tance of the K-point, when performing DFT calculations. 
The atomic forces in graphite are long ranged. Therefore, 
when using the finite-differences approach and DF T^^i^^i^^ 
only the phonons commensurate with the supercell are cal- 
culated correctly. In linear-response calculations, on the other 
hand, the implementation of the iiT-point can be more easily 
achieved, ^^'^^ however, the dynamical matrix at the i^T-point 
should be explicitly calculated and not be simply interpolated. 



TABLE I: Symmetry relations between space group and molecular 
notation for the space group F63 /mmc and the point group Deh. 
The corresponding coordinate system for C2V is shown in Fig. [5] 
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Phonon dispersion relations are often predicted from ab- 
initio DFT or from empirical force-constants (FC) calcula- 
tions. Empirical force-constants models in graphite have so 
far included up to 4th-nearest neighbors, in order to reproduce 
the overbending of the optical branch near the F point^^^ 
In the case of graphite, besides the details regarding the fre- 
quency values, both methods show differences in the shape of 
phonon branches, e.g., the position of the crossing of the ZO 
and TA modes near the M point. While in FC calculations 
a crossing of the ZO and TA branches between the F and M 
points is predicted, it is found in ab-initio results to take place 
between K and M. This probably stems from a fit to the only 
available TA mode from EELS experiments in Ref. [sl as our 
force-constants fit will show later (e.g. see Refl27b. 

Further differences between force-constants and ab-initio 
DFT calculations are found regarding the LA and LO 
branches near the M point: In DFT results, the LO-derived 



A. Fitting procedure 

To obtain the optimal fit of the set of force constants to 
the experimental data, we applied a variable neighbor search 
(VNS) method of global optimization^^, using the simplex 
method^^ for the local optimization subroutine. Basically, 
this is the least squares pr ocedure minimizing the average 
deviation A(f) = j^^J^i 1^^^^ - ^i(f)p between the 
^exp ^ experimental values, cj^^^, and the correspond- 
ing calculated frequencies, cj^ (f ), obtained by calculation with 
the trial values of the force constants f = (/i, . . . , /f). Due 
to symmetry, for each level of the 3 or 6 neighbor atoms the 
same triple of force constants can be used. We used stretch- 
ing, out-of-plane and in-plane force constants fi for each rel- 
evant pair of atoms Here, the direction of stretching cor- 
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FIG. 6: (Color online) Optical eigenmodes of graphene at the F point 
from the force-constants calculation. The experimental frequency 
values are given in brackets; they are taken from the data recorded 
closest to the F point. 



responds to the line that connects the center atom with the 
atom of appropriate level. The in-plane and out-of-plane di- 
rections are perpendicalur to this line and lie in the graphene 
layer or perpendicular to it, respectively. After transforming 
the stretching and in-plane force constants to a global basis, 
one obtains the dynamical matrix. The level of the relevant 
neighbors has been gradually increased until the satisfactory 
agreement (A < 0.23meV, with the greatest difference of 
Icj-^'P - cJi(f ) I ^8 meV for the LA branch in the F - re- 
gion, nearby the iiT-point) has been eventually achieved with 
included neighbors of up to the fifth level. The fifth level con- 
tains 24 neighbors of each atom, as there are 3, 6, 3, 6, and 6 
symmetrically positioned first to fifth neighbors, respectively, 
resulting in fifteen independent variational parameters fi . 



B. Force-constants results 

The optimized values of the force constants parameters are 
presented in Table HIl In Fig. [2] we show the phonon disper- 
sion obtained from these force constants in comparison to the 
experimental data. The largest deviations between the calcula- 
tion and the experiment occur for the optical phonon branches 
near the K point. This is probably due to the strong inter- 
action of the near-i^^-point phonons with electrons near the 
Fermi level, which is not included in a force-constants model. 
Including more than 4th-nearest neighbors of atoms (i.e., fif- 
teen independent parameters), however, gives a fairly good 
description of the local minimum of the TO-derived branch at 
the K point. 

Moreover, although FC calculations including only fourth- 
nearest neighbors provide a considerably good average fit to 
the experimental data, they lead to permuted frequencies of 
the LO and LA-derived phonons at the M point (M^ and 
M^), and to a crossing of the LA and LO branches within 
the - M region (K^ - M+ and - M^). Therefore, at 
least fifth-nearest neighbors are required for a good empirical 
description of the graphite phonon dispersion. 

The eigenvectors of all optical phonons from our force- 
constants calculation at the high-symmetry points F, and 
M are drawn in Figs. [51 [71 and[8l respectively. They are in 
agreement with calculations from a molecular-based approach 




FIG. 7: (Color online) Eigenmodes of graphene at the K point. For 
the degenerate modes, only one choice per energy is given. Phonon 
energies given in brackets are the IXS experimental values. For the 
symmetry notation see Table [J 



TABLE II: Force constants parameters for graphene, obtained from 
a fit to the experimental data, in eV/A^. 



Neighbor level 


Stretching 


Out of plane 


In plane 


1 


25.880 


6.183 


8.420 


2 


4.037 


-0.492 


-3.044 


3 


-3.016 


0.516 


3.948 


4 


0.564 


-0.521 


0.129 


5 


1.035 


0.110 


0.166 



of Ref. 28. For the degenerate modes, we show only one 
choice per energy; the remaining eigenvectors can be obtained 
by the symmetry-group projectors. 

Often in literature the molecular notation for the symmetry 
groups is used. Therefore, Table [D shows the relation between 
the spacegroup notation of P63 / mmc and the molecular nota- 
tion at the high-symmetry points F, and M, and the lines 
F - K - M (T) and F - M (S). The eigenvectors will help to 
choose the sample orientation in future IXS experiments. The 
scattering cross section is zero, if the direction of the atomic 
displacements and the momentum transfer in the scattering 
process enclose an angle of 90° . 

The above force constants parameters can in principle be 
used to calculate the phonon dispersion of carbon nanotubes, 
in particular for chiral nanotubes with a large number of atoms 
in the unit cell that require large computational effort, when 
calculated with first-principle methods. The existence of a 
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FIG. 8: (Color online) Eigenmodes of graphene at the M point (from 
lower to higher frequencies). Phonon energies given in brackets are 
the IXS experimental values. For the symmetry notation see Table U 

fourth acoustic mode (pure rotation of the tube about its axis) 
and the finite frequency of the radial breathing mode have to 
be taken into account, see Refs. i31ii32^33a . In addition, the 
different bond angles and lengths between the carbon atoms 
on the cylinder surface depending on the chiral index must 
be included. We expect that such an approach based on 



the empirical force constants of graphite will give an overall 
good description of the phonon bands of carbon nanotubes. 
In metallic nanotubes, however, the coupling of the F-point 
and i^T-point phonons to the electronic system will lead to 
different results for those modes, see for instance the Kohn 
anomalies and the frequency drop of the LO phonon in metal- 
lic nanotubes Therefore, where the phonon dispersion 
is modified due to strong interactions between the phonons 
and the electrons, a force-constants model might only give 
an emipirical description of the phonon bands, but should be 
tested by DFT calculations which take electron-phonon cou- 
pling into account. 

V. SUMMARY 

In summary, we presented the full in-plane phonon dis- 
persion of graphite determined by inelastic x-ray scattering. 
The overall shape of the phonon bands confirms previous ab- 
initio DFT calculations, if special care is taken for the high- 
est optical phonons near the K point. We showed that by in- 
cluding fifth-nearest neighbors, the phonon bands can be well 
described within a force-constants model. Previous empiri- 
cal models predicted only parts of the phonon dispersion cor- 
rectly, since experimental data in the — M region had been 
missing. The new force constants will also improve the mod- 
els of the phonon dispersion in carbon nanotubes. 
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